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O . We study the global dependence of large x F2 nucleon structure function data, 

(D ■ 

Q \ with the aim of providing a perturbative-QCD (pQCD) based, quantitative analysis 

' of parton-hadron duality. As opposed to previous analyses at fixed x, we use a 

^Sj ' framework in fixed W"^ . We uncover a breakdown of the twist-4 approximation with 

> : 

. a renormalon type improvement at 0{1/Q'^) which, by affecting the initial evolution 

\ of parton distributions, will have consequences for pQCD analyses also at large x 
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and very large Q 



One of the key challenges in Quantum ChromoDynamics (QCD) today is to formulate 



QhI a connection between the description of the hard, or short-distance, scattering processes, 
' which can be calculated in terms of quark and gluon degrees of freedom using perturbative 
, methods, and the physical asymptotic states, i.e. the spectra of hadrons which are not 

' calculable within perturbative-QCD (pQCD) and are in principle only remotely related to 
parton dynamics Yet, a substantial number of observations indicate that QCD has 
manifestly a dual parton-hadron nature. Recent investigations of hadronic jets at e~^e~ 
0) ep 0, and pp [^] colliders show that several infrared-sensitive features of the inclusive 
hadron distributions can be reproduced by a pQCD description of the parton shower down 
to Qo ~ ^QCD, Qo being an effective cut-off defining the onset of the non-perturbative 
regime and Aqcd being the scale of QCD. (Notice, however, that detailed angular correlation 
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measurements have been found to be in disagreement with available pQCD estimates 
Similarly, deep inelastic scattering (DIS) experiments at very low Bjorken x indicate that 
pQCD can describe the nucleon structure function F2 down to ~ 1 GeV^, where a smooth 
and fast transition to the non-perturbative regime occurs, leading eventually to the expected 
behaviour in the real photon scattering limit {F2{x,Q'^) —>■ for 0). The question of 

the role and nature of non-perturbative corrections naturally emerges. These are known to 
affect the cross sections for e~^e~ jet fragmentation and are also expected to be present, 
although they have not been isolated yet, in the deep inelastic data at low and low x [0. 

In this Letter we address yet another set of experiments, namely inclusive ep and eD 
scattering in the resonance region, where a pQCD description seems to hold in spite of 
the low values of the invariant mass squared, W"^, produced in the final state. Here, the 
observation of Bloom-Gilman (BG) duality [|], or the similarity between the behaviour of 
the resonance contribution to the nucleon structure function and DIS, can be formulated 
theoretically, as the equivalence between the moments of the structure function in the low 
kinematical region dominated by resonances and in the DIS one, modulo perturbative 
corrections and expectedly small power corrections Furthermore, because of the large 
body of highly accurate data that has been recently made available |Ty], the data in the 
resonance region can now be fitted to a smooth curve which, once evolved according to 
pQCD, will coincide with the DIS data. 

There has been a growing theoretical interest in this intriguing phenomenon, where the 
quasi-exclusive nature, or low inelasticity of the scattering process would be expected to 
hinder a direct observation of the partonic structure of the target. On the contrary, the 
measured structure functions appear to behave, in average, like parton distributions. The 
present analysis is motivated by the expectation that a more precise understanding of the 
mechanisms behind BG duality might provide a handle on the type of hadronic configurations 
that are present in a "semi-hard" regime, before confinement settles in. It is therefore 
of paramount importance to perform a detailed study of both the logarithmic corrections 
and of the size and nature of the power corrections in the pQCD expansion, in order to 
ascertain whether the apparent weak dependence of the data is coincidental, an artifact 
of the particular region under study, or a cancellation of Higher Twist (HT) terms, possibly 
understandable within parton-hadron duality models. The analysis conducted here involves 
a number of steps similar to recent extractions of power corrections from inclusive data 
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111 , [T2| , |T3|, |T4|, namely the form 



Q 

is adopted, where F^'^'"'^'''^^^'" (x, Q^) is the twist-2 contribution, including kinematical 
power corrections from the target mass (TMC); the other terms in the formula are the 
dynamical power corrections, formally arising from higher order terms in the twist expan- 
sion. Both j?pQ^^-^'^'^^^ and H can be extracted from the data at large x, by taking care of 
aspects of pQCD evolution peculiar to this region, and of TMC. Uncertainties are introduced 
at each step. In this Letter, we address them one by one with the important addition that 
we provide for the first time in the literature, both an analysis of the scale dependence of 
the resonance region (W'^ < 4 GeV^) and a combined study of both the resonance and the 
large W"^ data. 

We start by studying the question of the limit of applicability of pQCD in the low W^, 
large x domain, but at above the 1 GeV^ scale, so still amenable, in principle, to a 
pQCD treatment. Previous analyses used moments of the structure functions obtained at 
fixed values of and directly related to the OPE, circumventing the problem of dealing with 
the complicated structure of the bound states by providing a natural averaging procedure 
in the Mellin conjugate n 0. However, we have shown in |15| that, at lower values of Q^, 
such moments of F2 are affected by elastic scattering in such a way as to render higher twist 
contributions impossible to extract or interpret in the standard OPE language. We choose, 
therefore, an alternative method using data in fixed bins. 

The data in the resonance region |Ty] have been fitted to a smooth curve for F2^^\^), 



where ^ = 2x/{l+\/l + AM'^x'^ /Q"^) [jT6| is the Nachtmann scaling variable. The fit is applied 
separately to bins in invariant mass, W^, centered at: (1) = 1.6 GeV^, (2) = 2.3 
GeV^, (3) = 2.8 GeV^, (4) = 3.4 GeV^, respectively. The xVd.f. for these fits 
varies from 0.8 to 1.1. The uncertainty for the scaling curves is estimated to be better than 
10%, taking into account uncertainties in the experimental data, and in the averaging and 
fitting procedures. The curves and uncertainties in each W"^ bin are represented by the 
hatched areas in Fig. 1. Notice that the spectra at fixed require that vary for each 
spectrum as an increasing function of x, since = Q'^{x) = {W^ — M^)x/{1 — x). 

We utilize a quantitative comparison of the resonance region F2 curves thus obtained 
to the pQCD prediction for F2 at the same kinematics to extract potential higher twist 
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contributions to the structure functions. We assume that only valence quarks contribute to 
Fl^"\ at low < lOGeV^ and x > 0.3. PQCD evolution to = Q^{x) or at fixed W^, 
for the Non-Singlet (NS) distributions at next-to-leading-order (NLO), is given by: 



Jq'' 



Q'.2 271 



X f'^PgJ-,as{Q'''))qt\y,Q''% (2) 
J X y w J 



where q\ ^ = q^ — = q^, i = u,d. The expressions for the splitting function Pgg {z, as{Q'^))-, 



and for the corresponding endpoints at NLO in MS scheme can be found in [0, |T^; as is 
the strong coupling constant. The structure function, F2, is then obtained by convoluting 
(^) with the quark coefficient function, -B^"^ |l7|. We fix the values of the initial parton 
distribution functions (PDFs), at Q'^ ^ 0.4 — 1 GeV^, to the ones taken from NLO global 
fits to world data [^, and we solve the evolution equations directly in x space, with 

0:5 (M|) = .117. Notice in this respect, that the shape of the initial NS PDFs is practically 
constrained [^, at variance with the singlet and gluon distributions at low Q^, whose 
shape is strongly correlated with the value of as- The structure of perturbative evolution, 
including both coefficient functions and splitting functions for the NS part, can now be 
evaluated up to NNLO. Detailed studies of the impact of NNLO corrections and beyond, on 
the determination of power corrections for the NS structure functions, have been performed 
in ITT], O]: The question of whether these can "mimick" the contributions of higher twists, 



including the uncertainties due to the well known scale/scheme dependence of calculations, 
within the current precision of data is under intense investigation. Here, we single out the 
contributions that are expected to dominate the higher order perturbative predicitions at 
large x, namely powers of ln(l — z) terms, where z = x/y. These terms can be resummed 
to all orders. We perform the resummation in x space by replacing the scale with a 
z-dependent one, W"^ = Q'^{1 — z)/z p2[. It is well known that such a procedure introduces 
in principle an ambiguity in the evaluation of the running coupling at low values of the scale 
W'^, i.e. as z ^ 1 — A^/Q'^ (||2^ and references therein). This ambiguity is lessened in our 
analysis because at fixed W"^, A'^/Q^{x) is very close to zero. Our results for NS evolution 



and shown elsewhere |2J], are in good agreement with more recent resummation calculations 
performed in n space and anti-Mellin transformed as in p5| . 

Finally, we take into account target mass corrections (TMC). We use the expression in 
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0, obtained by a formal inversion of the Nachtmann moments [|T6]. Ambiguities in this 
procedure are expected to arise because of the neglect of higher order and higher twist 
corrections. To safely disregard these, we use a minimal criterion that only the kinematical 
points yielding low values of the parameter x^M^/Q^ in the TMC expansion are kept. 

Comparisons of pQCD+TMC predictions with the resonance average data are shown in 
Figures 1 and 2. A few comments are in order. In Fig. 1: i) TMC (full lines) modify 
substantially the pQCD behavior (dotted lines) rendering a better agreement with the data; 
a) the curve calculated using NLO pQCD at = 200 GeV^ shown for comparison in Fig. 
1(b) demonstrates the large effect of pQCD corrections above x ~ 0.2. 

In Fig. 2, we show the low data obtained from the fit of Jefferson Lab (JLab) and 
SLAC data in the resonance region, along with larger data from p6| , |27[| . Note that 
the data in the resonance region (circles) smoothly blend to the deep inelastic (stars) - 
another manifestation of BG duality. The curves correspond to our calculations including 
pQCD+TMC at NLO (dashes), and pQCD+TMC with resummation (full). The dots in each 
curve represent regions where TMC are uncertain. The effect we find is qualitatively similar 
to what found in 1^, in that over the range 0.45 < x < 0.85, higher order perturbative 
contributions, in this case large x resummation, improve the agreement with the data. 
Substantial discrepancies remain which we assume to be largely due to the dynamical, HT 
corrections to the structure function. If H{x, Q"^) is modeled similarly to previous extractions 



n|, 0, ig, one has: 

H{x, Q') = F|«^^+^^^^(x, Q')ChAx), (3) 

Eq.(^ is motivated by the lack of knowledge of the anomalous dimensions of the twist-4 
operators, a reasonable assumption within the precision of the data (see also |2^). Our 
analysis at fixed explained in the first part of the paper, enables us to extract Cht from 
the resonance region and from the large (DIS) region, separately. 

In Fig. 3(a) we show the coefficient Cht, Eq-(iD; extracted from: i) DIS data with W"^ > 4 
GeV^; ii) The resonance region, W"^ < 4 GeV^; in) Averaged over the entire range of W^. 



The figure also shows the range of extractions previous to the current one |jT3|, We notice 
that in all three cases, our values for Cht are smaller than the ones in [|13], because of 
the effect of large x resummation. We have checked that our results without resummation 
are consistent with a previous extraction using moments of the structure function . Most 
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importantly, while the large data track a curve that is consistent with the l/W"^ behavior 
expected from most models the low data yield a much smaller value for Cht and 
they show a bend-over of the slope vs. x, already predictable from a similar bend-over in 
the slopes at low W'^ in Fig. 2. 

Notice that this surprising effect is not a consequence of the interplay of higher order 
corrections and the HT terms, but just of the extension of our detailed pQCD analysis to 
the large x, low W'^ kinematical region. In fact, resummation produces an overall reduction 
of Chti that does not modify our results qualitatively. In order to ascertain whether the 
discrepancy between the low W"^ and large values of Cht are due to 0{1/Q'^) terms 
in the twist expansion, Eq.(0), which could become more important at low W"^, we have 
extracted for each resonance the quantity AH{x,Q'^), defined as 



ppQCD+TMC + Q2 



Aif(x,g2), (4) 



2 

where Cht{x) coincides with the value fitted at large W'^. From Fig. 3(b) one sees that 
Aif(x,(5^) is negative for all lower W"^ (< 3.4 GeV^) bins, as expected if a cancellation 
among higher order inverse powers were to occur, consistently with the requirement of 
parton-hadron duality. However, we uncover a non trivial dependence of this term: 
one can see a sharp change in the behavior of the higher mass resonances and the data in 
the N — A transition region which show a distinctively steeper fall with Q^. Furthermore, 
the high mass resonances seem to be in fair agreement with a simple fit inspired 

by Infrared Renormalon (IRR) calculations |^, which would produce a constant line at a 
decreasing height for increasing W"^ . The — A transition region shows a departure from 
this behavior which cannot be accounted for by fitting the first few terms of an inverse power 



expansion. A similar conclusion was reached in our analysis in moment space [|T^ . 

To summarise, the size of the power corrections obtained from the data at 0.45 < x < 
0.85, with W"^ > 3.4 GeV^, is comparable to the values obtained in recent analyses of DIS, 
with W"^ > 10 GeV^ [11, |l^, |l^, |2^, |2^, the inclusion of large x resummation producing 
a reduction of the coefficient Cht- At lower invariant mass, 1.9 < W"^ < 3.4 GeV^, an 
extra term of order 0{1/Q'^), with a negative coefficient is necessary to fit the data, in line 
with the asymptotic nature of the twist expansion and with current predictions from IRR 
calculations [3C]. A similar expansion does not reproduce, however, the scale dependence of 
data at even lower masses, 1.2 < W"^ < 1.9 GeV^, yet for values where pQCD would be 
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expected to apply. Therefore, the experimental observation of a rather flat dependence 
cannot be taken as a signature of parton-hadron duality for this region, but on the contrary, 
it shows the limits of applicability of this idea. 

Having this determination in hand, one can now attempt to provide theoretical interpre- 
tations. In particular, the fact that the pattern of dynamical power corrections to resonance 
production is comparable to the one for the large deep inelastic data, i.e. that power cor- 
rections remain small in the resonance region, suggests that color confinement is more likely 
to happen locally, with a smooth transition between partonic and hadronic configurations, a 
mechanism supported also by recent studies of hadron spectra in jet measurements [Q. This 
mechanism seems to break down in the N — A transition region, providing a threshold where 
cooperative effects from many partons dominate the structure function. More information 
on the locality of color confinement may be obtained by studying semi-inclusive processes 



where a more direct relation to cluster hadronization models can be established 24 
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FIG. 1: Comparison of NLO pQCD calculations (dashed lines), and NLO pQCD+TMC (solid 
lines) with the data on (hatched areas) at fixed values of W'^ = W^, vs. x: (a) W'^ = 1.6, 
(b) Ty2 ^ 2.3, (c) = 2.8, (d) = 3.4 GeV^. The dotted curve shows for comparison the 
DIS calculation, obtained at = 200 GeV^. The data are averaged with the procedure described 
in the text and reference. The pQCD curves were obtained using the GRV parametrizations for 
the NS distributions at the input value of = 0.4 GeV^ |1£], indicated by the arrows. Other 
parametrizations give similar results, starting from their input value of ~ 1 GeV^. Below 
= 0.4 GeV^, the similarity between the pQCD curves and the data in the resonance region can 
no longer be tracked by pQCD. 
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FIG. 2: Comparison of pQCD+TMC calculations at NLO (dashed lines) and with resummation 
(full lines), with current large x data. The full dots are in the resonance region, 1.3 < W'^ < 3.4 
GeV^; the open triangles correspond to W'^ < 1.3 GeV^. The dotted lines represent the regions 
where TMC contributions are uncertain. 
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FIG. 3: (a) Coefficient Cht, Eq.®, extracted from DIS data with > A GeV^ (full dots), from 
the resonance region, W'^ < 4 GeV^ (stars) and averaged over the entire range of W'^ (open dots). 
The shaded area summarizes extractions previous to the current one. A dotted line at zero is 
added to guide the eye; (b) AH, Eq.(^, extracted at fixed values of W"^ as described in the text, 
and plotted vs. Q^. The figure further elucidates a breakdown of the twist expansion at low W'^, 
already visible in (a). 



